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Abstract. In a recent study |arXiv:1011.3254j the contact process with a modified 
creation rate at a single site was shown to exhibit a non-universal scaling behavior with 
exponents varying with the creation rate at the special site. In the present work we 
argue that the survival probability decays according to a stretched exponential rather 
than a power law, explaining previous observations. 



1. Introduction 

In non-equilibrium statistical physics directed percolation (DP) is known to be the 
most important universality class of phase transitions into absorbing states, playing a 
role comparable to the Ising model in equilibrium statistical mechanics [IHZ]- One of 
the most intriguing properties of this class is its robustness with respect to the choice 
of the dynamical rules. Non-DP critical behavior is obtained only if the dynamical 
rules are changed substantially, e.g. by introducing additional symmetries, quenched 
randomness, or algebraically distributed long-range interactions. 

Recently Miller and Shnerb [8] studied the one-dimensional contact process with a 
single defect. The contact process is a version of DP with random-sequential updates 
following the reaction-diffusion scheme A 2A, A — and is often used as a simple 
model for epidemic spreading. Depending on the rate A for particle creation, it displays 
a second-order phase transition between a fluctuating phase and the absorbing state 
(empty lattice) taking place at the critical point Ac — 3.2978. Miller and Shnerb 
implemented a defect at a particular site, where the birth rate A is modified by the 
amount AA. Starting the process with a single active seed at this particular site, they 
observed that the survival probability Ps(t) measured at A = Ac exhibits a power law 
Ps{t) ~ t-^ with a varying exponent S depending on AA (cf. Fig. [1]). This is surprising 
since the variation of a rate in a single site can be considered as a boundary effect which 
usually does not modify the universality class of the process in the bulk. One should 




Figure 1. Left: Directed bond percolation with a defect. In the bulk each bond is 
open with probability p, while for bonds originating at the defect site (bold line) this 
probability is modified by Ap. The figure shows a particular configuration, where a 
cluster of active sites (blue dots) generated by an active seed at the defect site survives 
until time tg. Those bonds for which - if cut - the cluster would not have survived 
until to are called red bonds and are marked in red color. Right: Decay of the survival 
probability in a critical bond DP with a defect for various values of Ap. 



not confuse this with boundary induced phase transitions into an absorbing state [S], 
where the special dynamical rules in the boundary site generates the phase transition. 

In this paper we explain these findings. Considering the density of so-called red 
bonds at the defect site it turns out that the actual decay of the survival probability at 
criticality does not obey a power law, instead it is dominated by a stretched exponential. 
Our findings are confirmed by improved numerical simulations. 



2. Weak defect approximation 

In this section we derive an expression for the survival probability in DP with a weak 
defect, using the notion of so-called red bonds. Since red bonds are usually discussed in 
the context of directed bond percolation, we will use this realization of DP instead of 
the contact process throughout the paper. 



2.1. Red bonds 



Let us first consider DP without a defect. For a given cluster of sites connected by 
bonds, a bond is called 'red' if its removal would break up the cluster into two pieces. 
In the context of the present problem, where we study the survival of a DP cluster 
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generated from a single seed, the red bonds of a cluster surviving until time t are those 
which, if removed, would turn the cluster into a non-surviving one (see Fig. [T]). 

In DP, the density of red bonds is closely related to the off-critical behavior of 
the model. To see this let us consider an ordinary bond-DP process without defects 
at criticality p = Pc and ask the question how the survival probability will respond to 
an infinitesimal variation of the global control parameter p. Suppose that for a given 
realization of open and closed bonds a seed at the origin has generated a cluster that 
survives until time t. Let us then slightly decrease the percolation probability by Ap 
and ask for the probability that the cluster turns into a non- surviving one. Clearly, by 
decreasing p (using the same set of random numbers) some of the conducting bonds of 
the cluster will be removed. If the removed bond happens to be a red bond, the cluster 
breaks up into to pieces and becomes non-surviving. Therefore, the response of the 
survival probability to a variation of p is directly related to the actual number of red 
bonds. 

More specifically, let Ntot{t) be the total number of red bonds of the cluster and 
consider a small change of the percolation probability p — > p+Ap by Ap < in a system 
without defect. Then the survival probability will respond to lowest order by 

P,{t,p, + Ap) = Ps{t,pc){l + aAp)^-'W + . . . (1) 

with some model-dependent constant a, so that to leading order the change is given by 

P,(t,Pe + Ap) ~ t~^ exp(a ApiVtot(t)). (2) 

Since a change of the percolation probability will lead to a breakdown of the algebraic 
behavior of the survival probability at a time scale of the order |Ap|~'^ii, we can 
immediately conclude that the total number of red bonds in a cluster scales as 
A^tot(t)~t+'/''". 

2.2. Red bonds in the center 

Let us now turn to directed percolation with a single defect. Here the situation is 
analogous to the one described before, the difference being that p is not changed globally 
but only at the site from where the cluster was generated. Therefore, the response of 
the survival probability will be determined to leading order by the number of red bonds 
originating from this particular site, meaning that one simply has to replace the total 
number of bonds Ntotit) by the number of red bonds in the center, denoted as Nc{t). 
Since we can expect the density of red bonds to be roughly constant in the interior of the 
cluster, the (time-integrated) number of red bonds in the center will be approximately 
proportional to A'tot(^) divided by the typical spatial volume of the cluster V ~ t*^/^. 
This means that Nc{t) is expected to scale as 
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Figure 2. Numerical verification of formula The graph shows the same data 
as in Fig. 1 multiplied by i*"" , where dcS is the effective survival exponent averaged 
in the time interval from 10^ to 10^ time steps. In this representation, clean power 
laws would appear as horizontal lines. However, the numerical data (black) exhibits 
a clear curvature depending on the strength of the defect. The suggested stretched 
exponential ([5|) fitted to the data by adjusting the constants ci and C2 captures this 
curvature convincingly. 



Using this scaling law we expect the survival probability to behave to leading order as 
Ps(t) ~ exp(a A|>ti/"ii-'^/") (4) 
with some non-universal constant a. 



3. Numerical results 

3.1. Survival probability 

The argument presented in the previous section implies that the survival probability in 
1+1-dimensional DP with a single defect should decay as 

P,(t)^Cir°-1^95g^p^^^ ^-0.05586)^ (5) 

provided that Ap is small and negative. In order to verify this formula, we performed 
extensive list-based seed simulations in a virtually infinite one-dimensional system, 
averaging over 10^ to 10^ runs depending on Ap and over 10^ time steps. Note that the 
simulations in Ref. [5] cover only 30.000 time steps, - a much smaller range where the 
curvature of the data in loglog representation was probably not easy to see. 
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Cl 


C2 


C2/Ap 


-0.256 


0.230 


0.078(5) 


2.12(5) 


-8.3 


-0.128 


0.195 


0.29(2) 


1.05(3) 


-8.2 


-0.064 


0.177 


0.52(1) 


0.53(2) 


-8.3 


-0.032 


0.168 


0.70(1) 


0.26(2) 


-8.1 





0.159 


0.95(1) 


0.01(2) 




+0.032 


0.150 


1.24(2) 


-0.27(2) 


-8.4 


+0.064 


0.142 


1.63(3) 


-0.53(3) 


-8.3 


+0.128 


0.126 


2.68(5) 


-1.01(4) 


-7.9 


+0.256 


0.094 


6.8(1) 


-1.94(5) 


-7.6 



Table 1. Numerical estimates of the effective exponent of the survival probability 
(average slope of the data without error margin) and the fit parameters ci and C2 for 
various values of the defect strength Ap (see text). 

To demonstrate the non-algebraic behavior, we first determined the effective 
exponent (5eff by an ordinary power-law fit over four decades of time ranging from 10^ 
to 10^ Monte Carlo steps (see Table 1). Then we multiplied the data by f^^^: In a 
loglog representation a clean power law would give a horizontal line so that even tiny 
deviations are much easier to see. In fact, as shown in Fig. |2l there is a systematic 
curvature depending on the strength of the defect. 

A non-linear fit of Eq. (|5]) over four decades reproduces this curvature accurately. 
It is surprising that the formula, which was derived in the limit of small and negative 
Ap, works very well for positive Ap and even up to |Ap| = 0.256, which is more than 
30% away from the critical point. The corresponding estimates of the constants ci and 
C2 are listed in Table 1. 

As an additional check, we verify the prediction that the second fit parameter C2 
should be proportional to Ap. In fact, the ratio C2/Ap listed in Table 1 is approximately 
constant except for large values of Ap, where we see a small deviation. 

3. 2. Number of active sites 

Besides the survival probability, seed simulations provide data for the total number of 
active sites N{t) averaged over all runs and the mean square spreading from the origin 
R^{t) averaged over all surviving tuns. In critical DP without defects, these quantities 
scale as N{t) ~ and B?{t) ~ t^/^, where rj = \/z — 25 and z = i^\\/i^±- 

As shown in Fig. [31 the quantity N{t) shows similar deviations like Ps{t) if the defect 
strength Ap is varied. However, the quotient N{t)/ Ps{t), which can be interpreted as 
the number of active sites averaged over surviving runs only, seems to exhibit clean 
power laws. This quotient divided by the expected power law t^~^^ is shown on the right 
hand side of the figure. As can be seen, small deviations of next leading order are still 
visible. Nevertheless the quotient N{t)/Ps{t) approximates a power law much better 
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Figure 3. The number of active sites N{t), the survival probabihty Ps{t), and the 
ratio N{t)/ Ps{t) in seed simulations of DP with a single defect at the origin. 

than N{t) and P,(t) in itself. 

3.3. Density at the defect site 

DP with a defect can also studied by simulations with homogeneous initial states. For 
example, we may ask how the density at the defect site decays as a function of time 
starting with a fully occupied lattice. It turns out that this density decays in the 
same way as the survival probability analyzed above. This is no surprise since rapidity 
reversal, a symmetry implying Ps{t) and the density at the "seed site" starting with fully 
occupied lattice to decay with the same exponent, holds also in presence of a defect. A 
proof of this relationship is given in Appendix A. 

4. Summary 

In this work we have shown that a single defect in l+l-dimensional DP does not change 
the universality class of the process, opposed to the conclusion obtained in Ref. [8]. 
For weak defects with Ap/p ^ 1 we derived a formula for the survival probability by 
considering the density of red bonds at the defect site. 

We note that in the extreme case Ap = —p the defect site acts as an absorbing 
boundary. Even in this case the universality class of the bulk process does not change, 
instead one obtains an interesting boundary critical behavior with a new surface critical 
exponent p^lTTj. Away from this point, this boundary critical phenomenon breaks 
down and the defect site has no influence in the renormalization group sense, instead it 
manifests itself as a correction of certain quantities such as the survival probability in 
form of a stretched exponential. 
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Appendix A. Time-reversal symmetry 

It is well known that directed bond percolation without defects has the special property 
that the time reversal symmetry is exact. It can be checked easily that this proof can be 
generalized to percolation probabilities with any type of quenched disorder, including 
the single defect special case. 

In order to make sure that this symmetry is not an artifact of directed bond 
percolation but valid in any realization of DP, we show here that the time-reversal 
symmetry also holds in the contact process for transition rates A depending on the 
lattice position. The calculation below is similar to the proof of time-reversal symmetry 
for the homogeneous contact-process presented in [T^ . 

As in Ref. [12] we define the creation and annihilation operators 

a'i\si) = Si\si - I) , a||si) = (1 - + 1), (A.l) 
where Sj = 0, 1. The operator H for the CP is given by 

H = $^[(al - l)ai] + - l)ala]a,]. (A.2) 

i <ij> 

The density of particles at some site i is given by 

p,(t) = {0\a,Vmt)) (A.3) 

where V = Ylji^ + Oj)- Note that only configurations with a particle at i "survive" the 
action of the operator ciiV. If the initial condition is a fully occupied lattice we have 

p,(t) = {0\a,Vexp{-Ht)l[a]\0). (A.4) 

j 

The pseudo-Hermitian operator H = VHV~^ obeys the relation [12] 

W = aHa, (A.5) 

where o" = (j^ = a"^ = Yljl^j^^]] = Ylji^ ~ ^j)- Note that a possible 

inhomogeneity of the creation rate has no influence on the last equality. 

Finally, Eq. flA.4p implies that 

p,{t) = {0\aiexp{-Ht)V^O) = {0\Vaexp{-Ht)aa^^\0). (A.6) 

Since ^'''ay^lO) = |0) and o"a||0) = — a||0), we obtain the final result 

p,(t) = -(0| exp{-Ht)V"'al\0) = 1 - (0| exp(-^t)al |0) = P,(t). (A.7) 

This proves the time reversal symmetry in the contact process with position-dependent 
creation rates. 
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The fact that time reversal symmetry holds provides a direct way of observing that 
a single defect site does not alter the universality class. If the density of active sites 
(averaged over the whole lattice) starting with a fully occupied lattice is measured we 
expect a clean DP decay exponent, which imply in the single defect site generating a 
transient effect that is felt near enough the boundary. 
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